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Abstract 

We consider the contribution to K ^ 27i decays from the non- 
diagonal s ^ d quark transition amplitude. First, we calculate 
the most important part of the s ^ d transition, the so-called 
self-penguin amplitude ~ G_f«s, including the heavy top-quark 
case. Second, we calculate the matrix element of the s ^ d tran- 
sition for the physical K 27t process. This part of the analysis 
is performed within the Chiral Quark Model where quarks are 
coupled to the pseudoscalar mesons. 

The CP-conserving self-penguin contribution to — 27r is 
found to be neghgible. The obtained contribution to e'/e is sen- 
sitive to the values of the quark condensate < qq > and the 
constituent quark mass M. For reasonable values of these quan- 
tities we find that the self-penguin contribution to e'/e is 10-15 % 
of the gluonic penguin contribution and has the same sign. Given 
the large cancellation between gluonic and electroweak penguin 
contributions, this means that our contribution is of the same 
order of magnitude as e'/e itself. 
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1. Introduction 



In general, non-leptonic AS = 1 decays are described by an effective 
Lagrangian at quark level ||, || 

CiAS = l)=J2C^Q^ . (1) 
i 

In such effective Lagrangians, the heavy mass scales are integrated out and 
their effects are contained in the Wilson coefficients Cj. The quark opera- 
tors Qi involve the three light quarks q = u,d, s. The coefficients will in 
general include short distance QCD effects calculated by means of perturba- 
tion theory and the renormalization group equations (RGE). In our notation 
the Wilson coefficients Cj also contain Fermis coupling constant Gp and the 
relevant Kobayashi-Maskawa factors Xq = Vsq^dq CP-conserving (g = u) 
and CP-violating (g = t) AS" = 1 transitions respectively. 

Within this standard procedure, one usually omits operators containing 
{i'j ■ D — nig) , by appealing to the QCD equations of motion for quark fields 

{t^ -0-711^)^0, (2) 

where is the covariant derivative containing the gluon field, and ruq the 
appropriate quark mass. This procedure corresponds to going on-shell with 
external quarks in quark operators. Certainly, quarks are not on-shell in 
hadrons, especially not in the (would-be Goldstone) mesons vr and K. 

In a series of papers 0, ||, it has been pointed out that the effect 
of such operators are not zero in general, and that they in particular will 
contribute to processes like — >■ 27r , 27 , and B — *■ 27. In order to be 
consistent, such genuinely off-shell effects should go to zero in the limit when 
some typical bound state parameter goes to zero, which is indeed found to be 
the case. In our case this parameter is the constituent quark mass of the light 
u, d, s quarks. Numerically, the non-zero off-shell effects are non-negligible 
in some cases. 

In this paper we will study a quantity which by construction vanishes on 
the s- and d- quark mass-shells after regularization, namely the non-diagonal 
s ^ d self-energy transition ||^, 0, ||, [1^]. For off-shell quarks, however. 



this renormalized non-diagonal self-energy is still non-zero and could a priori 
be relevant for physical amplitudes, e.g K 2tt. The unrenormalized s —>■ d 
transition corresponds to an effective Lagrangian 



C'^^ = d{ai-f ■ DL + bL + cR)s , (3) 

where a, b, c are divergent quantities given by the loop integrations, and L, R 
are the left- and right-handed projectors in Dirac space. In the Standard 
Model Lagrangian there are no direct s ^ d transitions. Consequently, one 
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demands that the renormahzed self-energy vanishes on the respective mass 
shells of the s- and d-quarks, by adding the necessary counterterms. It should 
be emphasized that this defines the physical s and d quarks in the presence of 
weak interactions. The renormalized self-energy corresponds to an effective 



Lagrangian of the form ITT 



C^^ = -Ad{n ■ D - md){h ■ DR + MrR + MLL){i-f D - ms)s . (4) 

In momentum space, A is a finite, slowly varying (logarithmic) function of 
the VT-boson and heavy quark masses and the external s- and d- quark 
momentum (k) squared. To be explicit, A is obtained from an expansion 
of a in (^: a = + k'^A + ■ ■ ■. The quantity Oq, which is divergent and 
independent of fc^, is removed by the counterterms. The quantities Ml^h 
depend on the masses of the s- and cZ-quarks. 

In the limit rris^d 0, which we will work, the Lagrangian is 

Ci{ms,d ^0) = -A d{i^ ■ DfLs . (5) 

In the pure electroweak case, there is a strong suppression of the CP- 
conserving amplitude owing to the GIM- mechanism |12|, and one finds 



that A is of order Gf^I/ M^r. In the CP- violating case this strong GIM 
suppression is relaxed due to a top-quark with a mass of the same order 
as the W-boson. Even more important, when perturbative QCD to lowest 
order is added, one obtains an unsuppressed charm-quark contribution ~ 
Gpotsilogm^'^ [0, H. A Feynman diagram for this contribution is shown in 
Fig. 1 and is named self-penguin because the gluon from a one loop penguin 
diagram is reabsorbed by the external s- or d-quarks. 

If one applies the perturbative QCD equations of motion as in (^, we 
observe that 0. However, we should emphasize that C^^ is obtained 

within perturbation theory. Therefore, it is not allowed to apply for 
quarks strongly bound in vr, K. An analogue within QED would be to replace 
(iby 

(7 ■ - m) ^ 67 ■ , (6) 

where the perturbative photons are contained in the covariant derivative, 
and A^ represents the binding Coulomb forces [|13|. Also, in QED the renor- 
malized electron self-energy is zero on the electron mass-shell, but still it 
gives an important contribution to Lamb-shift. As a gedanken experiment, if 
the s- and ci-quarks were bound by electromagnetic forces, the contribution 
to [ds) —>■ 27 from (|^) would be proportional to the binding energy of the 
system, that is, of the order aem/rB, where is the Bohr radius. Being 
small in QED, one expects off-shell effects to be bigger in strong interactions 
where bound state effects are more important. Thus, physical effects for 
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Figure 1: The self-penguin diagram. There is also a counterpart where the 
gluon is absorbed by the rf-quark. 



K ^ 2ti decays from C^^ could be obtained, and one should explore possible 
consequences for the A/ = 1/2 rule and for e'/e. 

To calculate the matrix elements of quark operators between physical 
hadronic states is in general a difficult task, and one normally uses various 
models or assumptions. In this paper we will use the Chiral Quark Model 
(xQAf), advocated by many authors |T3]. This model consists of the ordi- 
nary QCD Lagrangian and an extra term which induces meson-quark cou- 
plings that makes it possible to calculate matrix elements of quark operators 
in terms of quark loop diagrams. The model is thought to be applicable 
for quark momenta below some scale of the order A,^, the chiral symmetry 
breaking scale. We find this model very interesting because it reproduces the 
chiral Lagrangian terms for strong interactions to good accuracy to order p^. 
Furthermore, it has been shown to be useful in calculating weak hadronic 
matrix elements O, ITH . 



2. The s ^ d self-penguin transition 

The renormalized non-diagonal s d off-diagonal self-energy has the 
following form in agreement with (^) 



M{s d) = -idT^dsS = -ik^ d'j ■ kLsA{k) 



(7) 



where k is the four momentum of the external s- and d- quark. In the pure 
electroweak (EW) case we obtain for one fiavour q = u,c, t: 



EW 



2A„-^ 1 + 



167r2 



dx— In 

k^ 



1 - 



x(l — x)k^ 



(8) 
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where = xM^r + (1 — x)m^, My/ is the ly-boson mass and m = niq is the 
mass of the quark running in the loop (g = u,c,t), gw is the Ty-couphng, 
and \q is the appropriate KM factor. The term ~ rn? / (2M^) is due to 
unphysical Higgs exchange in the Feynman gauge. For small k"^ we obtain 

G 1 

Aew = -^^[A«(r„-rc) - At(rc-n)] , (9) 
where the quantities are given by 

for q = u,c, t. As already mentioned, the pure electroweak CP-conserving 
s ^ d transition amplitude is very suppressed due to the GIM mechanism. 
A priori one expects a relaxed GIM mechanism for the CP-violating case 
involving the heavy top quark. The top quark contribution is, in Feynman 
gauge, dominated by the unphysical Higgs contribution. As a curiocity, we 
find a "super-GIM" mechanism: It turns out that the top-quark contribution 
exactly cancels the charm contribution {{r^ — rt) —>■ 0) in the limit —>■ 
and — i> oo. Therefore the GIM-cancellation is somewhat stronger than 
expected in the pure electroweak CP-violating case. Numerically, we find 
(rc - n) = -0.167 - (-0.124) = -0.043 for mt = 180 GeV. 

As shown originally by Shabalin 0, the self-penguin diagram has only a 
weak, logarithmic GIM suppression to order agGp. To obtain the two loop 
self-penguin we need the penguin {s ^ d + gluon) loop for arbitrary external 
momenta. By direct calculation this is divergent, and before we insert it into 
the next loop to get the self-penguin, it has to be properly regularized. For 
this purpose we consider the Ward-identity 

g,r^, = S,,(fc + g)-S,,(fc) , (11) 

in order to determine the necessary counterterms. In our case F^^ corresponds 
to the penguin one loop (the strong coupling Qs not included). Writing 

r:^,(fc,g)=T-7.^, (12) 

the most general tensor compatible with the Ward identity (|rT|) can be writ- 
ten 



rpO 



{k + qfA{k + q)- k''A{k)\ ^ + {k + qfA{k + q) g'' 



+ Boiq'-q'' - q'gn + B,[q'^k'' - q ■ kg^""] 
+52[fc'^fc^ - ^k'^q"] + B^te'^'-'^^k^qf, . (13) 
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Comparing this with the raw resuh for V^^, we find the regularized expres- 
sions: 



B, = 2gl. f dx f ' dyf^ [m^ - x{l - x)q^ - y{l - y)P + 2xyq ■ k] ' .(14) 

For W- exchange only, the f/s are given by /o = 2x{l — x); fi = y{2x — 
l);/2 = — 2|/2 and = —y, respectively. There are similar expressions 
for the unphysical Higgs case. We obtain the following expression for the 
renormalized self-penguin contribution: 

Mis-^d,g',) = ~2zglgl,dt^e J J^:^S{k + q)V'^,,{k,q) s , (15) 

where t"" is a colour matrix. To obtain the result compatible with (^, we 
expand the d- (or s-) quark propagator S{k + q) and FJ^^ in powers of the 
external momentum k and keep in total the order (/c)^ only. Then we obtain 
(for q = c,t): 

M{s ^ d,g',) = -^3^X^ekW'^^,Ls (-^) V, , (16) 



3M2, •"■^ Vl67r2 



where the dimensionless quantity: 
JO Jo 



[M2 - x(l - x)A2 



yx ^ ?/2(3x — 1) x{l — x) 



M2 (1 - x)M2 3M2 

»1 rl-x 



+3m^ dx dy 
Jo Jo 



[(-3/2 + 2x)y + x{l - x)]/n(^(^) 



[M2 - a;(l - a;)A2] 
yix + 3y) 3y\l-y) ^ 2y^ xjl^x) ] ^^^^ 



M2 s(l-x)M2 (l-a;)2M2 3M2 

includes the unphysical Higgs contribution ~ 3m2. Since we have omitted 
any current mass for the s- and rf-quark and expanded to third order in the 
external momentum k, one has to introduce a lower cut-off scale A. In order 
to match the low-energy description in terms of the xQM, its magnitude will 
be taken as the chiral symmetry breaking scale A-^ = 0.83 GeV. 

It should be noted that in the previous calculation 0], only the case 
m2 <^ was considered. In this case only the Bq term contributes to 
order Gp, while contributions from the other i?j's are suppressed by m? JM'^ 
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with respect to the Bo contribution. Including the heavy top quark, all 
contributions have to be taken into account. 

Combining the equations and (|T^, we find the self-penguin contri- 
bution to the quantity A in (^: 

= ^ ^ f [A.(F. - F,) - A,(F, - F,)] , (18) 

where he = TritH") /N^ = {N^ - l)/{2Nc) is a colour factor, and A^^ is the 
number of colours. Eqs. (p!5|)-([T^) is the contribution from the diagram in 
Fig. 1 (including the unphysical Higgs contribution). Taking into account the 
counterpart where the gluon is absorbed by the d-quark, one gain a factor 
two which is included in ([T8|) . Numerically, we find for m = rric = 1.4 GeV 
and A = 0.83 GeV: 

Fc~49.0 . (19) 
Similarly, for m = rrit = 180 GeV and A = 0.83 GeV, we obtain 

Ft ~ -10.7 . (20) 

For the c-quark case one obtains the analytical result within the leading 
logarithmic approximation 

= -M^W - -IH^)? + (---)ln(^) - -ln(^) . (21) 
The leading part of this result was obtained in [§] for the charm quark case, 



while the general result ( |T^ ) is new. In ( pl| ) the term ~ —4/6 is coming from 
—B^q^g"^ in ([T3|) and the rest from B^q^q'^. As shown in ||l^ |18|, leading 
logarithmic QCD corrections can be taken into account by folding the well 
known Wilson coefficients C± H into an integral over virtual momenta p. 
Then our result takes the form 

(22) 

where p{p^) = [C+(p^) +C_(p^)]as(p^)/(27r). Because the operator in (|^) has 
zero anomalous dimension, there is no extra factor attached to p(p^). The 
linear combination (C+ + C_)/2 is the colour octet combination of the co- 
efficients which occurs because the gluon within the self-penguin is a colour 
octet. The C± 's are given by ratios of at different scales to a power 
determined by the anomalous dimension of the corresponding four quark op- 
erators Q±. Thus the integral in ( p2D can be performed analytically, and the 
result obtained in terms of at the scales A,mc,mb, and Mw respectively. 
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Numerically, we find Fc — 3.25. For the top quark contribution, we expect 
essentially no QCD corrections except for the running of as, and we take 

Ft = ^^^Ft - -0.37 , (23) 

71 

approximately the same as the pure EW contribution. The total A to be 
used in (||) is Asp + Aew- (There are other contributions of order agGp, but 
with a very small numerical coefficient. Note also that Aew is dominated 
by loop momenta of order Mw where C± ~ 1, and there will be no further 
QCD corrections in the leading logarithmic approximation). 



3. The Chiral Quark Model 

In the Chiral Quark Model (xQM) [Q, chiral-symmetry breaking is 
taken into account by adding an extra term to ordinary QCD: 

C^ = -M{qZ^^qR + <m^qL) , (24) 



where q = {u, d, s). The constant M in (|2^) is interpreted as the constituent 
quark mass, expected to be of order 200-300 MeV. The quantity S contains 
the Goldstone-octet fields Tia- 

= exp{iJ2^aTTa/ f) , (25) 

a 

where the A^'s are the Gell-Mann matrices, and f = f^w = MeV is the 
pion decay constant. The term £^ introduces a meson-quark coupling. This 
means that the quarks can be integrated out and the coefficients of the various 



terms in the chiral Lagrangian are calculable 



The xQM should be applied for momenta lower than the scale of chiral- 



symmetry breaking, = 2nfTryQ/Nc, where A*";, = 3 is the number of 
colours. Divergent integrals are regularized, for instance with dimensional 
regularization [0. An alternative is to use some ultraviolet cut-off Ajj of 
the order A^. Because /t^, entering the meson- quark coupling ~ M'j^/ f, is 
also given by a quark loop diagram for vr W^(virtual), one has: 

(0) _ A^ ~ TT^ a 2 

- +mjr^^n J' ^^^^ 

where < ^G"^ > is the two gluon condensate, and the dots indicate higher 
gluon condensates. In the end both / and f^^^ will, in the limit m„ — > 0, 
be identified by but at intermediate stages one might need to distinguish 
them for technical reasons. The dimensionless quantity /t^ has the leading 
behaviour ~ ln(A^/M^) in a cut-off type regularization and ~ r(e) ~ 1/e 
within dimensional regularization {D = 4 — 2e being the dimension of space). 
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Figure 2: Typical diagram for non-diagonal s — > c? self-energy induced K — »• 
27T transition. 

The quark condensate can be written as 

NcM ^ 1 tts ^9 
<?,>=^C,-^<-G^>+-. (27) 

The quantity Cq depends on the regularization prescription. Using a cut- 
off A, one finds Cq = —A? + ln(A^/M^) + ■ ■ ■ , and within dimensional 
regularization, C„ = — M^r(— 1 + e) + ■ ■ ■. 



4. The K ^ 271 amplitude from s ^ d. 

Within the xQ^i the K ^ 27t amplitude obtained from the s d 
transition in is determined by diagrams like the one in Fig. 2. We find 
the following result: 



TT TT 



V2 



s = -77^(^1 - ml) AM' 



< qq > 
M 



+ VI 1 



3M2 
2A^ 



(28) 

we have also calculated 
TT transition amplitudes 
These calculations show 



where A is the quantity in (J^). For completeness 
the virtual K vacuum (which is zero) and K — 
in addition to the physical K ^ 2tt amplitude, 
that in (|) gives an octet contribution that can be simply added to the 
Qq contribution up to OijP'). To order 0{j)'^) we expect to obtain different 
contributions from and C^g due to the extra derivatives in 

We want to compare the result (|2^) with the corresponding amplitude 
for the penguin operator Qq : 



M 



X J 



(29) 
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We observe that the ratio between the amphtudes in (|28| ) and (^) goes to 
zero in the hmit where the constituent quark mass M ^ 0, as expected. 

Numerically, the ratio between the self-penguin and the Qq contribution 
is less than 5 % in the CP-conserving case. For e'/e, the situation is different 
because there is a large cancellation between the gluonic and electroweak 



penguin contributions |]T9|, |2y, |I6|. In the operator language, it is a large 
cancellation between the Qq and the Qg contributions. We have compared the 
s ^ d self-energy contribution with the standard contribution in the CP- 
violating case. Using, in the notation of Ci = —{Gp/ V^){^uZi — XtV 
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with 1/6 = —0.137 at the renormalization scale /i = A = A^; M ~ 220 MeV 
and < qq >i/3~ —220 MeV [0, we find that in the CP- violating case, the 
s — >■ contribution is 10 to 15 % of the Qq contribution and has the same 
sign. This means that after the partial cancellation between the Qe and the 
Qs contributions has been taken into account, the s —>■ d contribution is of 
the order e'/e itself. 

5. Discussion 

We have considered the s —>■ d self-energy contribution to ^ 27r within 
the chiral quark model. Within this model matrix elements between the light 
mesons (vr, K) of quark operators may be calculated in terms of quark loops. 
Therefore this model provides a means to calculate effects of off-shell quarks. 

In the literature, terms which vanish by using the perturbative equations 
of motion are omitted from the analysis from the very beginning P, ^, |]. 
But to neglect such effects would, within QED, be analogous to discarding 
Lamb-shift. Moreover, such effects are larger in QCD where the bound state 
effects are stronger. Still, also in strong interactions such effects are relatively 
small, being proportional to M^/A^ ^ 0.07 and vanishing in the free quark 
limit M — > 0. They might, however, become non- negligible in cases where 
there are cancellations between the potentially largest terms, as shown in 
this paper for e'/e , and for i^' ^ 27 p where the potentially large pole 
contributions cancel when the Gell-Mann-Okubo mass formula is used. 

Numerically, we have found a small effect for the CP-conserving K 271 
amplitude. However, for the CP-violating K ^ 271 amplitude and thereby 
for e'/ e the effect is non-negligible, being of the order e'/ e itself, because of the 
strong cancellation between the Qg and Qs contributions. It would of course 
have been an advantage if were included in the basis of a complete next 
to leading QCD corrections pOI, but this would have required a three loop 



calculation. The self-penguin contribution drives e'/e more to the positive 
side. However, due to the uncertain value of the quark condensate it is still 
difficult to make a very precise prediction of e'/e within the chiral quark 
model |T6[]. 
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